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1.  INTRODUCTION 


The  Rydberg  states  of  the  molecule  have  long  been  of 
interest.^  As  a  result  of  the  axial  symmetry  of  the  molecule^  in  contrast 
to  the  spherical  symmetry  of  atoms,  one  expects  a  wide  variety  of  Rydberg 
series  for  N^  •  Because  of  the  complexity  of  the  N2  spectrum, 
some  predicted  Rydberg  series  have  not  been  observed,  and  there  are 
many  observed  lines  as  yet  unidentified. 

13  + 

In  this  paper  we  shall  present  theoretical  calculations  of  the  ’  Z  , 
13+13  13 

’  ’  n  ,  and  ’  Rydberg  states,  all  converging  to  the  ground 

2  +  + 

X  Z  state  of  the  N„  ion.  Numerous  theoretical  treatments  of  N-  Rydberg 

g  2  2  ^  “ 

2-5 

states  have  appeared  in  the  literature.  All  these  works  attempted 

to  construct  an  effective  one-electron  Hamiltonian  to  solve  for  the  Rydberg 

2 

orbitals  as  is  done  in  this  work.  In  this  regard  Duncan  and  Damiani 
approximated  the  core  orbitals  by  one-center  functions  (centered  at  the 

midpoint  of  the  two  nuclei)  to  facilitate  the  calculation.  Betts  and 

3  5  3 

McKoy  used  a  model  potential  in  the  Hamiltonian.  They  used  as  basis 

functions  a  set  of  Caussian-type  orbitals  (GTO)  centered  at  various  points 

4 

along  the  nuclear  axis.  Lefebvre-Brion  and  Moser  constructed  the 


Hamiltonian  by  using  the  occupied  orbitals.  Their  basis  functions  for  the 
Rydberg  states  are  the  Slater-type  orbitals  centered  at  the  two  nuclear  sites. 

In  many  ways  our  approach  parallels  that  of  Ref.  4.  However,  we  take 
the  view  that  a  Rydberg  state  must  reflect  a  physical  situation  in  which  an 
electron  moves  in  a  nearly  central  field.  Thus,  we  begin  by  treating  a 
quasi-spherlcal  molecular  system  (centered  at  the  midpoint  of  the  molecule) 
in  which  a  Rydberg  electron  is  characterized  by  atomic-llke  orbitals  with 
angular  momentum  (t,m).  Since  this  Idealization  is  not  fulfilled  for  the 
real  molecular  Rydberg  states,  we  make  allowance  for  t-mixlng  to  obtain 


■ 


3 

the  final  results.  Formal  theory  on  the  Rydberg  states  Is  extensively 

6  6,7 

discussed  by  Mulllken.  A  detailed  analysis  of  the  molecular  spectra 

requires  consideration  of  such  interactions  as  the  coupling  between  the 

electronic  motion  with  molecular  rotation.  However,  we  restrict  our 

scope  to  the  electronic  part  of  the  problem,  and  calculate  the  wave 

functions  at  the  equilibrium  separation  R  *  1.1  X  of  the  X  state  of 


11.  METHOD  OF  COMPUTATION 


2  +  + 

The  Rydberg  states  converging  to  the  ground  X  state  of  N2  ion 
have  the  electronic  configuration, 

+  2  + 

Since  the  core  (N2)  is  in  a  state,  Rydberg  orbitals  of  the  type 

ka  ,  ko  ,  kr  ,  and  kn  give  rise  to  the  ,  and  ^*^11 

g‘  u’  g’  u  *  g'  u»  g*  u 

states  respectively.  The  corresponding  wave  functions  are  constructed  as 
antisymmetrized  products  of  these  one-electron  molecular  orbitals  (MO) 


eg*. 


s  s 

+  |lo  (rj^)a  lo  (r2)6...3o  (rjj_2)6  (rjj)a|], 


where  |...|  represents  a  determinant,  and  a,B  are  the  two  spin  functions. 
The  Hamiltonian  of  the  system  is  (in  atomic  units) 

where  r^,  are  the  two  nuclear  sites  and  z(=7)  is  the  nuclear  charge. 

We  adopt  the  frozen-core  approximation,  namely,  we  assume  that  the  N2  core 

orbitals  are  not  affected  by  the  addition  of  the  Rydberg  electron.  Hence 

the  la  ,  2o  ,  lo  ,  2a  ,  In  ,  and  3o  core  orbitals  in  Eq.  (2)  are  obtained 
g  g  u’  u*  u’  g  M  ^  / 

8  2  ^ 

by  the  self -consistent-field  (SCF)  procedure  for  the  X  Z  state  of  N_. 

g  2 

9 

In  this  SCF  calculation  the  GTO  are  used  as  the  basis  functions,  and  the 
core  orbitals  are  expanded  by  GTO  centered  at  the  two  nuclear  sites. 

As  for  the  Rydberg  orbitals  we  express  them  in  terms  of  basis  functions 
which  are  centered  at  the  midpoint  of  N„,  viz.. 


TBI  JWPJiP.'ff.ipj  jijy.»r.'i.t^iiij;>.  |ij|,.j  I 
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■'Ry 


’ 


with 


(4) 


(5) 


The  notation  is  standard  except  that  n  is  replaced  by  k.  Our  use  of  one-center 
basis  functions  for  the  outer  electron  differs  from  the  works  of  Refs.  3  and 
4  in  which  the  basis  functions  consisted  of  molecular  orbitals  constructed 
from  prechosen  atomic  orbitals  judged  to  be  suitable  by  the  respective 
authors. 


A.  Basis  Functions 

Following  the  standard  notation  of  the  Hartree-Fock  theory,  w«;  designate 
the  electron-nucleus  and  electron-electron  interaction  by  V  and  the  electron 
exchange  operator  by  W.  For  a  diatomic  molecule  if  we  use  the  spherical 
harmonics  as  the  angular  part  of  the  basis  functions,  as  in  Eq.  (5),  the  axial 
nature  of  the  diatomic  molecule  makes  V  and  W  diagonal  in  m  but  not  in  i. 

This  is  reflected  in  the  i-summation  that  appears  in  the  expansion  of  the 
wave  functions  according  to  Eq.  (4) .  To  ensure  fast  convergence  of  the 
k-summatlon  in  that  expansion,  we  choose  as  the  functions  the 

solution  of  a  quasi-spherlcal  molecular  Rydberg  system.  In  other  words 


the  are  obtained  by  solving  the  integro-dif ferential  equation 


Ki+l) 


dr 


2  ■  ^k£in^ 


(6) 


where 


6 


''«..f.‘'>  ■  ’'to''>’'f.<'>  + 1^1'^ 


+2  /dr  Yj^*(r)Y^,j^(r)  Ej  fj  /  ♦*(r*)  jr-r' |”%j(r')dr* .  (7] 

The  suiranatlon  over  j  covers  all  the  core  orbitals  with  the  respective 
occupation  numbers  .  The  exchange  potential  Is  an  Integral  operator,  viz., 

‘mr’' 

>«  [{r'^'Vr'V,  (r’)dr*  +  r*"  /”r*'’''\, (r')dr’ } 

OK.x.iik  r  KJCiB 

xr'/d;’  Yj.g(;’)Yjj,^(r')^*(?')l.  (8] 

where,  again  the  summation  Is  over  the  core  orbitals  4>j  with  appropriate 
numerical  coefficients  f ^ .  The  detailed  numerical  procedures  of  computing 
Eqs.  (7)  and  (8)  are  shown  elsewhere. Here,  we  merely  note  that  Eq.  (7) 
is  computed  with  no  approximation,  but  the  summation  over  t  in  Eq.  (8)  is 
an  infinite  series,  which  is  truncated  after  five  terms. 

Except  for  the  exchange  term,  we  follow  very  closely  the  procedure  of 
Herman  and  Sklllman^^  to  compute  numerically  A®  exchange 

term  we  do  not  use  the  Slater  approximation,  but  adopt  a  scheme  which  has 


been  applied  successfully  to  electron-scattering  problems. 


That  is. 


as  we  solve  Eq.  (6)  iteratively  for  etc.  for  the  n-th 

and  (n+l)st  sets  of  solutions,  we  approximate  the  integral  operator  as 

ypCn+l)  ^  ]/F^"^  }  , 


except  for  the  initial  cycle  where  we  set  WF' 


0*  The  function  WF' 


is  first  computed  by  using  the  solution  of  the  previous  cycle  f'‘  as 
shown  in  Eq.  (8),  then  divided  by  F^'*\  Thus,  Eq .  (6)  remains  a 


7 


homogeneous  equatlcxi  as  in  the  calculation  without  exchange  or  with  exchange 
treated  by  Slater's  approximation.  In  the  limit  of  convergence 
in  the  self-consistent  iteration,  Eq.  (9)  is  justified.  We  have  not 
encountered  any  serious  numerical  problems  associated  with  dividing 
by  The  mesh-sizes  used  are  shown  in  Table  I. 

In  this  manner,  we  computed  the  radial  basis  functions  F,  .  for  k  =  1-11. 
i  =  0-9,  and  m  *  0,1.  The  exchange  coefficients  f'  appropriate  for  the  singlet 

•J 

states  are  used.  No  separate  set  of  basis  functions  is  calculated  for  the 
triplet  states. 

In  Tables  Il-V  we  list  the  values  of  which  are  the  eigenvalues  of  the 

quasl-spherical  Rydberg  molecule.  For  large  values  of  k,  we  expect  to 

approximate  fairly  well  the  Rydberg  levels.  In  this  paper,  of  course,  we 
do  not  take  and  as  the  solutions  of  Rydberg  states.  Instead  we  use 

the  functions  as  basis  functions.  We  expect  these  functions  to  outperform 

most  other  arbitrarily  chosen  basis  functions. 

B.  Dlagonalization 

From  Eqs.  (2)-(5),  we  have  a  system  of  equations  for  the  coefficients 

with  the  matrix  elements 


\£m,k’£’m  =  F^^^(r)[  2  "  2^ 

dr  r 


(11) 


(12 


(r)y,,  (r)  drdt. 

%  B 


Because  we  used  the  frozen-core  approximation,  c  obtained  from  Eq.  (10)  does 
not  yield  wave  functions  that  are  orthogonal  to  the  occupied  orbitals.  To 
ensure  the  orthogonality  we  make  the  following  modification.  For  example, 
in  the  case  of  Z  states,  we  construct  a  new  set  of  orthogonal  basis  functions 
by  using  the  c's  obtained  from  Eq.  (10),  i.e., 

MW 


V  g 


k«  kl,v  ktm  11, 


(13 


where  i|i^(o  )  are  the  three  occupied  orbitals  (lo  ,  2o  ,  3o  ) ,  and  the  index 


g  g  g  g 

V  indicates  the  collating  sequence  of  eigenvalues.  The  x  functions  are  then 


used  as  basis  to  construct  the  Hamiltonian  and  overlap  matrices  H*  and  S'. 
Solution  of  new  secular  equation 


(H’  - 


t'S')  b 


0. 


gives  the  orthogonalized  Rydberg  functions,  e.g.. 


(14; 


♦/V  ■ 


E  b  X  (o  ) 

V  V,  g 


kt 


C,  a  U 

kl  p  iclm 


(r)  +  ).  p  =  1,2,, 


(is; 


We  have  solved  eight  different  sets  of  secular  equations  corresponding  to 

the  ,  and  ^’^11  states.  The  MO  <(>.  in  Eq,  (15)  are  drawn 

from  the  occupied  orbitals,  (la  ),  (2o  )  (3a  ),  (la  ),  (lx  )  according  to 

g  g  g  u  u 

the  symmetry  of  each  secular  equation.  The  eigenvalues  obtained  from 
Eq.  (14)  are  shown  in  Tables  VI  -  XIII.  Here,  as  in  Tables  II-V,  we  present 
six  digits  for  the  eigenvalues  in  order  to  see  how  closely  the  c' -eigenvalues 
approach  the  X-eigenvalues  and  the  hydrogenic  levels  with  increasing  f.. 


III.  RESULTS  AND  DISCUSSION 


A.  The  eigenvalues 


Inasmuch  as  the  A  eigenvalues  In  Eq.  (6)  for  large  k  may  be  taken 
as  the  first  approximation  to  the  energy  levels  of  the  Rydberg  states.  It 
is  instructive  to  analyze  the  values  of  X  in  terms  of  the  effective 
potentials  V  and  W  in  Eq.  (6).  The  index  k  specifies  the  number  of  nodes 
of  which  is  equal  to  k-i-1.  The  values  of  X  shown  in  Table  II-V 

reveal  the  following  features  in  comparison  with  the  energy  levels  of  the 
H  atom  which  are  Included  in  the  last  column  of  Table  IV: 

(1)  The  values  of  with  £  ^  4  are  very  close  to  the  corresponding 

hydrogenlc  levels,  suggesting  that  these  basis  functions  may 
already  be  a  fair  representation  of  the  respective  states. 

(2)  The  values  of  with  £  =  0,1,2  show  no  resemblance  to  the 

corresponding  hydrogenic  levels. 

(3)  For  £  =  3  the  X-values  of  the  (kf ) (£  =  3)  series  appear  to  be 
displaced  nearly  by  one  unit  in  k  from  the  hydrogenic  levels. 


whereas  the  members  of  the  (kf)  11^  series  are  close  to  the 
hydrogenic  levels  of  the  same  k. 


The  finite  separation  of  the  nuclei  is  an  obvious  cause  of  departure  from  the 
central  field.  Figs.  1-3  illustrate  the  difference 


AV,  „  =*  V„  „  -  (-  1/r).  (16) 

£m ,  £ni  £m ,  £m  ^ 

Because  of  the  strong  attraction  by  the  nuclei,  AV  shows  a  potential  cusp 

at  the  nuclear  site  (r  =  R/2  =  1.05  a  ).  AV  is  negative  for  r  >  0.3  a  . 

N  o  °  n 

It  covers  a  range  of  r  'v  0.6  -  1.5  a^  and  becomes  negligible  at 
r  -v  2.0  a^.  The  position  of  the  cusp  is  evidently  dictated  by  the 
geometry  of  the  inol ocular  ion,  but  the  depth  of  V  depends 


tv  *;  vs: 


I  -.Jj  ^  -Jf 


not  only  on  the  core  orbitals  but  also  on  the  syimoetry  (i.ml  olC  the  Rydberg 

state  under  consideration.  For  instance  o-type  Rydberg  orbitals  in 

general  overlap  more  favorably  with  the  charge  density  of  the  core 

than  do  u-type  Rydberg  orbitals.  This  is  reflected  by  the  larger  |aV| 

for  the  I  states  since  all  the  Rydberg  states  considered  in  this  paper 
+  2  + 

have  a  common  1^2  (  ^  )  core.  The  attractive  perturbation  AV  must  be 

2 

weighed  against  the  repulsive  centrifugal  term  il(t+l)/r  in  Eq.  (6), 

For  i  4  the  centrifugal  barrier  completely  overpowers  AV  and  keeps  the 

^ktm  away  from  the  origin  that  they  have  virtually  no 

overlap  with  AV.  The  corresponding  A’s  should  be  very  close  to  the  hydrogen 
2 

value  of  -’-ik  .  For  smaller  S,  the  perturbation  AV  shifts  the  nodal  points 
of  relative  to  the  corresponding  hydroge'nic  wave  functions  as 

discussed  in  Refs.  6  and  13.  Moreover  the  strongly  attractive  potential 
well  of  V  has  a  number  of  deep  bound  states  which  correspond  roughly  to 
the  core  orbitals  of  N2.  These  deep  bound  orbitals  add  further 
complexities  to  the  nodal  structure  of  the  Rydberg  orbitals  through  the 
orthogonality  requirement.  Hence  we  find  a  wide  range  of  shift  of  A, 


from  the  hydrogen  limit  -h'k-  for  £  _<  3. 

The  effect  of  electron  exchange  may  be  estimated  perturbatively  by 
using  the  hydrogen  radial  function  R  ;  namely,  we  have  computed 


“kw  ■  / 

For  the  (ks)^i;^  series  the  ratio  |w/AV|  ranges  from  0.049  (k=4)  to  0.018 

(k=ll);  for  (kd)^E''',  from  0.0047  (k=4)  to  0,0021  (k=ll);  for  (kp)^^'’’, 

8  u’ 


11 


from  0.15  (k*A)  to  0.061  (k“ll);for  (kd)"*^!!  ,  from  0.0063  (k-4)  to 
0.0029  (k-11) ;  and  for  (kp)^n^,  from  O.IA  (k»A)  to  0.052  (k-11) .  For  the 
higher-angular  momentum  states,  the  electron  exchange  Is  quite  negligible 
(10~^  or  less). 

B.  The  Rydberg  Levels 


The  energy  levels  of  the  Rydberg  states  are  obtained  from  Eq.  (14). 

Due  to  the  orthogonality  requirement,  the  energy  roots  associated  with  the 

deep  levels  of  the  effective  potential  are  eliminated.  Equation  (14) 

allows  various  (kj.)  members  to  mix  within  a  symmetry  of  A  (or  m  since  the 
+  2  -f 

core  is  in  a  state) .  Except  for  a  few  cases  there  is  always  a 

dominant  (ki)  member  in  each  Rydberg  state  wave  function.  Thus  we  label 

a  Rydberg  state  by  n  which  is  taken  as  the  k-value  of  the  dominant  member. 

We  also  retain  the  label  £  even  though  it  is  not  strictly  a  good  quantum 

number.  The  Rydberg  state  energies  are  summarized  in  Tables  VI  -  XIII. 

13+  13  + 

By  far  the  strongest  mixings  occur  in  the  (ns)  ’  E  and  (nd)  *  E  series 

R  S 

and  the  mixing  coefficients  for  the  singlet  states  are  shown  in  Tables 

XIV  and  XV.  Admixture  with  (k£)  members  of  £  ^  4  are  much  smaller  and 

omitted  from  these  tables.  For  the  (ns)^!^  series  the  identity  of  n  is 

not  clear  until  we  reach  n  *  7.  We  simply  assign  n  =  6,5,'**'  to  the 

lower  roots  of  the  series  consecutively.  The  next  group  of  series  shows 

an  appreciable  mixing  with  the  adjacent  k  values  but  no  mixing  with  the 

13+  13 

different  £  values.  This  group  includes  the  (np)  '  E^,  (nd)  *  JI^,  and 
1  3 

(np)  ’  states.  For  illustration  abbreviated  tables  are  presented  of 
the  singlet  states  of  this  group.  The  mixing  coefficients  with  different 
£*  =  £  ±  2  are  at  most  0.011,  and  usually  much  less.  The  rest  of  the  states 
(£  ^  4)  are  entirely  atomic  in  character  with  the  dominant  coefficients 
i  0.9999.  A  comparison  of  the  X-eigenvalues  in  Tables  11  -  V  and 
with  the  corresponding  e*-eigenvalues  in  Tables  VI  -  XIII  also  confirms  this. 


The  mixing  patterns  for  the  triplet  states  are  qualitatively  similar  to 
those  of  the  singlet  counterparts  and  are  not  given  here. 

We  now  Interpret  the  energy  levels  In  Tables  VI  -  XIII  in  terms  of 
the  quantum  defect  6 


2(n-6n,)‘ 

The  m-dependence  Is  Implied  through  the  symmetry  of  Rydberg  states.  Many 

theoretical  analyses  have  been  made  on  the  quantum  defect  of  atomic 
14 

Rydberg  states.  Characteristically  the  quantum  defects  of  an  atom 

depend  on  1,  but  are  nearly  constant  of  n  within  an  f-serles.  Further,  If 

*'max  largest  of  the  angular  momenta  of  the  core  orbitals,  the 

quantum  defect  is  very  small  for  those  Rydberg  states  of  1  >  1^^^ . 

Although  a  general  analogy  to  the  atomic  case  clearly  exists,  such  a  simple 

characterization  of  the  molecular  Rydberg  states  must  be  tempered  with 

caution.  In  the  present  problem  of  N^,  the  quantum  defect  6^  clearly 

shows  an  2. -dependence  as  In  the  atomic  case.  But,  In  contrast  to  the 

atomic  case,  varies  considerably  with  n.  The  largest  variation  of 

6  ,  with  n  is  found  with  the  (ns)^!^  series  in  which  6  .  varies  from  0.718 
nf  g  nl 

(n=4)  to  0.413  (n=ll);  (nd)^n  ,  0.871  (n=4)  to  0.641  (n=ll); 

g 

and  (np)^E^,  1.581  (n=4)  to  1.309  (n=ll) .  These  states  are  characterized 
by  strong  mixing  of  the  basis  functions  as  shown  In  Tables  XIV,  XVI  and 
XVII.  The  next  group  shows  a  moderate  variation,  which  includes 
(np)^Il  for  which  fi  ranges  from  0.660  (n=3)  to  0.596  (n=ll) ;  and  (nd) 

U  TIJL  g  * 


1.0449  (n=4) 


to  1.0325  (n=ll);  and  (nf)^Ej^,  1.00678  (n=5)  to  1.00916 


(n=ll) .  These  states  have  relatively  small  mixing  among  the  basis  functions. 
In  the  case  of  the  (nd)^E^  series  an  appreciable  mixing  is  found  between  the 


basis  functions  of  the  (k,d)  and  (k-l,s)  series  (see  Table  XV),  but 


.  py  mj  VK w.' wirupuiru  w,  n  -V 


energetically  those  two  groups  of  basis  functions  are  very  similar.  Other 

states  of  large  I  show  markedly  smaller  quantum  defects.  This  Is  In 

accord  with  the  concept  of  penetratlng/nonpenetratlng  atomic  orbitals.^ 

[The  core  orbitals  may  be  decomposed  into  the  constituent  l-components . 

From  this  analysis,  the  charge  densities  are  found  to  be  mainly  confined 

tol£3:  lo  (42.6%),  2a  (92.5%),  3o  (94.8%),  lo  (60.0%)  2a  (96.9%) 
g  g  g  u  u 

In  (98.0%).]  However,  even  when  the  quantum  defect  is  small.  6  .  Is 

still  not  entirely  constant  within  a  series.  A  very  similar  situation 

exists  with  the  triplet  states  studied  in  this  paper.  Tliijs,  we  are  not 

able  to  support  the  qualitative  predlction^^  that  6^^  is  constant  of  n 

for  the  molecule.  Since  the  exact  solutions  of  are  available, 

Mulllken^^  analyzed  the  quantum  defect  of  as  a  function  of  the  internuclear 

separation.  At  R  =  2.0  a^  a  comparable  variation  of  with  n  is  found 

in  H^.  This  is  consistent  with  the  present  results.  Moreover,  the 

R-dependence  of  6^  may  be  understood  in  terms  of  the  combined  effect 

of  the  angular  momentum  barrier  and  the  potential  cusp  4V(r)  at  r^j  =  R/2. 

Unfortunately,  because  of  this  complexity  a  generalization  of  the  quantitative 

features  is  difficult  for  molecules. 

The  (nd)^L^  and  (nf)^Z^  series  have  a  quantum  defect  slightly  larger 

than  one.  If  we  shift  the  n  scale  by  one  unit,  these  two  series  would 

appear  to  have  very  small  quantum  defects.  The  same  shift  would  make 

the  apparent  quantum  defects  of  (np)^I^  smaller  and  those  of  (nd)^n 

^  g 

negative.  Care  should  be  taken  to  allow  for  this  possible  ambiguity  in 
the  assignment  of  n  when  comparing  theoretical  Rydberg  levels  with  experiments 
The  Rydberg  states  of  are  important  to  the  optical  properties  of 


the  atmosphere  as  transitions  among  many  of  these  states  are  in  the 
infrared  region.  Application  of  the  selection  rules  At  =»  ±1,  AA  =  0 
.and  tl,  and  AS  =  0  to  the  Rydberg  levels  reported  in  T.ablcs  VI  -  XIII 


lA 

yields  a  vast  number  of  transitions  of  wavelength  In  the  range  of  4-20 
ym  which  are  displayed  In  Figs.  4  and  5.  All  transitions  are  drawn  to  the 
same  height  and  no  consideration  was  given  to  their  relative  intensities. 

Since  we  have  considered  only  levels  through  n  *  11,  Figs.  4  and  5  do  not 
represent  a  complete  listing  of  all  possible  transitions  in  the  4-20  vnn  region 

To  illustrate  the  richness  of  the  Rydberg  spectrum  we  mark  in  these 
figures  the  wavelengths  of  the  limiting  hydrogenic  transitions  of  the 

Rydberg  levels  involved.  Figures  4  and  5  are  based  on  consideration  of 

electronic  states  with  no  reference  to  molecular  vibration  and  rotation. 

The  vibrational  potential  curves  of  all  the  Rydberg  states  considered 

here  have  virtually  the  same  shape  as  the  potential  curve  of  their  parent 
+  2  + 

^2  ^  ^g^  state  and  therefore  have  identical  vibrational  frequency. 

Furthermore  the  vibrational-level  selection  rule  for  electronic  transitions 
between  such  Rydberg  states  is  Av  =  0  according  to  the  Franck-Condon 
principles.  Thus  the  vibrational  levels  have  no  effect  on  the  transition 
frequencies.  The  rotational  levels  do  produce  fine  structures  in  the 
electronic  transitions  replacing  each  line  in  Figs.  4  and  5  by  a  band. 


IV.  COMPARISON  WIIN  OTHER  WORKS 
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The  energy  levels  of  those  Rydberg  states,  which  may  be  compared  with 
other  works,  are  compiled  in  Table  XIX.  In  accordance  with  the  previous 
works^  the  energy  levels  are  converted  to  show  the  term  values  above  the 
N-(X^i:'*^)  ground  state  based  on  the  fact  that  N^(X^J:^)  is  15.581  eV  above 

^8  2  g 

N  (xV). 

2  g 

4 

The  theoretical  calculations  of  Lefebvre-Brion  and  Moser  are  similar 

to  the  present  work.  Both  calculations  are  quite  rigorous  within  the 

framework  of  the  frozen-core  approximation.  In  describing  the  Rydberg 

orbitals,  however,  they  approached  it  from  the  separate-atom  picture, 

whereas  we  use  the  united-atom  picture  as  the  starting  point.  Therefore, 

the  reliability  of  the  present  work  increases  with  the  increasing  quantum 

numbers.  The  model  potential  employed  in  the  calculation  of  Ref,  3  does 

not  reflect  the  details  of  the  potential  of  the  one-electron  Hamiltonian 

discussed  in  the  previous  section.  Although  the  model-potential  approach 

affords  a  great  deal  of  computational  simplicity,  it  does  not  Include  the 

2  +  2  2  + 

provision  of  distinguishing  the  ionic  states  (X  ,  A  11^,  B  E^)  to 

which  the  Rydberg  states  converge,  or  the  spin  multiplicity  (singlet  or 

triplet)  of  the  Rydberg  states.  Nevertheless,  all  three  sets  of  theoretical 

calculations  agree  reasonably  well. 

The  experimental  data  are  summarized  in  the  extensive  review  article 

by  Lofthus  and  Krupenie.^  In  some  cases  the  experimental  values  cited 

17  18 

in  Table  XIX  are  taken  from  the  original  papers  ’  as  annotated.  A 

general  agreement  is  seen  between  the  theory  and  experiment.  Some 
19 

refinements  of  the  experimental  work  have  been  published  recently.  The 
new  data  agree  with  the  earlier  ones  to  within  typically  five  digits.  This 
is  below  the  estimated  uncertainty  of  our  theoretical  values. 


Most  of  the  observed  Rydberg  levels  are  confined  to  the  (np)^!^ 
and  (np)^n^  series.  Since  the  3a^  orbital  has  52Z  (l.“0)  and  43%  (1*2) 
components.  It  would  seem  that  (nf)^!^  and  (nf)^n^  states  may  appear  In 
the  absorption  spectra  analogously  to  the  (np)^r^  and  (np)^n^  states. 
However,  no  experimental  evidence  has  been  found  so  far  as  to  their 


existence  or  identity.  In  this  connection  we  may  mention  the  recent 


identification  of  the  k  II  state  as  a  Rydberg  state  converging  to  the 

R 

2  +  +  X8 

X  (N^) .  The  term  value  14.100  eV  is  close  to  14.191  eV  of  the  present 


calculation. 


17 


V.  SUMMARY  AND  CONCLUSIONS 

The  spectrum  of  the  N2  molecule  is  indeed  extremely  complex.  In  this 
paper  we  have  dealt  with  the  electronic  part  of  the  problem.  We  have 
placed  the  emphasis  on  the  calculation  of  high  Rydberg  states,  therefore, 
the  united-atom  approach  Is  a  natural  choice.  For  high-angular-moiDontum 
states  the  molecular  core  may  be  treated  as  a  quasl-spherlcal  one.  For 
low-angular  momentum  states,  mixing  of  the  eigenstates  of  the  limiting 
quasi-spherlcal  system  is  necessary.  The  agreement  with  other  theoretical 
calculations  and  experiment  is  reasonable  and  satisfactory. 

The  finite  separation  of  the  nuclei  adds  a  new  dimension  to  the 

problem  that  is  not  present  in  an  atomic  case.  It  appears  that  the 

properties  of  Rydberg  series  may  vary  greatly  from  one  molecule  to  another 

and  from  one  series  to  another  within  the  same  molecule.  A  rough  estimate 

of  the  electron-exchange  effect  is  about  15Z  of  the  Coul  mb  effect  or 

less.  In  the  frozen-core  approximation  adopted  here,  we  take  the  core 

+  2  -f 

orbitals  as  those  of  the  N,  (X  Z  )  ion  and  neglect  any  orbital  readjustment 

^  g 

due  to  the  addition  of  the  Rydberg  electron.  Such  a  charge  redistribution 
is  often  referred  to  as  the  core  polarization  which  is  sometimes  treated 
by  Introducing  the  polarization  potential  -a/2r^  where  a  is  the  polarizability. 
For  the  case  of  atoms.  Van  Vleck  and  Whitelaw^®  have  shown  that  the  use  of 
the  polarization  potential  is  valid  only  under  several  restrictive  conditions. 
To  fully  address  this  point  and  other  pertinent  questions  such  as  electron- 
correlation  energy  would  require  computations  on  the  level  of  configuration 
interaction.  Although  the  multiconfiguration  self-consistent-field  method 
has  been  applied  to  the  ground  and  low  excited  states  of  molecules,  it  would 
be  an  enormous  undertaking  to  incorporate  the  multiconfiguration  scheme  to 
the  Rydberg  states.  At  this  time  a  more  practical  direction  to  improve 


the  present  work  Is  to  allow  the  core  orbitals  to  relax  under  the  Influence 
of  the  Rydberg  electron.  The  core-orbital  relaxation  should  be  more 
Important  for  the  lower  Rydberg  states.  This  is  consistent  with  our 
observation  that  the  lower  Rydberg  states  show  larger  discrepancy  between 
theory  and  experiment  in  Table  XIX.  Furthermore  inclusion  of  core-orbital 
relaxation  by  a  variational  procedure  would  lower  the  calculated  Rydberg 
levels  making  a  change  in  the  right  direction  toward  the  experimental  values 
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Table  1.  Mesh  size  of  Integration  in  a 


0.0 

0.0125 

1.5 

1.5 

0.0250 

2.5 

2.5 

0.0500 

5.5 

5.5 

0.1000 

12.5 

12.5 

0.2000 

356.5 

and  are  the  starting  and  terminating  points  respectively  with  the 
step-size  6r. 
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1  +  a 

Table  11.  X-eigenvalues  of  E  (m  -  0)  type  basis  functions  . 

g 


2.92183(0) 

1.59652(-1) 

6.55096(-2) 

3.53493(-2) 

2.20659(-2) 

1.50711(-2) 

1.09426(-2) 

-2.21840(0) 

-5.60753(-2) 

-3.l4l38(-2) 

-2.00714(-2) 

-1.39264(-2) 

-2.00297(-2) 

-1.39061(-2) 

-1.02149(-2) 

-l,02075(-2) 

■8.30418(-3) 

-1.02263(-2) 

-7.8l968(-3) 

-7.81472(-3) 

■6.51646(-3) 

-7.82672(-3) 

-6.17785(-3) 

-6.17435(-3) 

-6.17345(-3) 

■5.24949(-3) 

-6.18248(-3) 

-5.00361(-3) 

-5.00107(-3) 

-5.00039(-3) 

■4.31904(-3) 

-5.00684(-3) 

-4.13492(-3) 

-4.13300(-3) 

-4.13249(-3) 
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X  ^  fl 

Table  III.  X-eigenvalues  of  (m  ■  0)  type  basis  functions  , 


k 

t  -  1 

£  -  3 

£  -  5 

£  “  7 

£  -  9 

2 

-3.53198(0) 

3 

-1.95374(-1) 

4 

-6.52276(-2) 

-2.25576(-l) 

5 

-3.46552(-2) 

-3.13678(-2) 

6 

-2.16090(-2) 

-2.00570(-2) 

-l,38980(-2) 

7 

-1.47799(-2) 

-1.39207(-2) 

-1.02098(-2) 

8 

-1.07502(-2) 

-1.02236(-2) 

-7.81628(-3) 

-7.81390(-3) 

9 

-8.17186(-3) 

-7.82535(-3) 

-6.17545(-3) 

-6.17377(-3) 

10 

-6.42207(-3) 

-6.18174(-3) 

-5.00188(-3) 

-5. 00063 (-3) 

-5.00027(-3) 

11 

-5.17998(-3) 

-5.00639(-3) 

-4.13359(-3) 

-4.13268(-3) 

-4.13239(-3) 

®See  the  footnotes  in  Table  II. 
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Table  IV.  X-eigenvalues  of 


n  (m  1)  type  basis  functions 

8 


k 

i  -  2 

i  -  4 

£  -  6 

£  -  8 

Hydrogenic** 

3 

-6.7A709(-2) 

-5.55556(-2) 

4 

-3.84349(-2) 

-3,12500(-2) 

5 

-2,39647(-2) 

-2.00253(-2) 

-2.00000(-2) 

6 

-1.62267(-2) 

-1.39034(-2) 

-l,38889(-2) 

7 

-1.16817(-2) 

-1.02132(-2) 

-1.02072(-2) 

-1.02041(-2) 

8 

-8.80125(-3) 

-7.81856(-3) 

-7.81455(-3) 

-7,81250(-3) 

9 

-6.86535(-3) 

-6.17706(-3) 

-6.17424(-3) 

-6.17340(-3) 

-6.17284(-3) 

10 

-5.50324(-3) 

-5. 00304 (-3) 

-5.00098(-3) 

-5. 00037 (-3) 

-5.00000(-3) 

11 

-4.50910(-3) 

-4.13448(-3) 

-4.13293(-3) 

-4.13248(-3) 

-4.13223(-3) 

*See  the  footnotes  in  Table  11. 

b  2 

Energy  levels  of  H  atom  (=  -  l/2k  ) 


>  .V.N  .%.%■%  ,v 


3 


5 


2 

-6.63839(-l) 

3 

-8.58743(-2) 

4 

-4.22486(-2) 

-3.13557(-2) 

5 

-2.52510(-2) 

-2.00572(-2) 

6 

-l,68054(-2) 

-1.39231(-2) 

-1.38971(-2) 

7 

-1.19918(-2) 

-1.02260(-2) 

-1.02092(-2) 

8 

-8.98732(-3) 

-7.82732(-3) 

-7.81590(-3) 

-7.81383(-3) 

9 

-6. 98617 (-3) 

-6.18330(-3) 

-6.17518(-3) 

-6.17372(-3) 

10 

-5.58631(-3) 

-5.00766(-3) 

-5.00166(-3) 

-5.00059(-3) 

11 

-4.56879(-3) 

-4.13796(-3) 

-4.13346(-3) 

-4.13265(-3) 

-5.00024(-3) 

-4.13237(-3) 


3 

-1.06273(-1)*’ 

4 

-4.64080(-2) 

-5.72547(-2) 

5 

-2.64274(-2) 

-3.19593(-2) 

-2.00298(-2) 

6 

-1.71486(-2) 

-2.03465(-2) 

-1.39062(-2) 

7 

-1.20559(-2) 

-1.40793(-2) 

-1.02150(-2) 

-1.02078(-2) 

8 

-8.94882(-3) 

-1.03180(-2) 

-7.81987(-3) 

-7.8l498(-3) 

9 

-6.90857(-3) 

-7. 88498 (-3) 

-6.17811(-3) 

-6.17465(-3) 

10 

-5.49274(-3) 

-6. 22097 (-3) 

-5.00386(-3) 

-5.00134(-3) 

11 

-4.46081(-3) 

-5.03262(-3) 

-4.13517(-3) 

-4.13328(-3) 

-6.17376(-3) 
-5. 00067 (-3) 
-4.13278(-3) 


11 
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Table  VII. 


1  a 

Energy  levels  of  £•  type  Rydberg  states 


n 

.-1^ 

t  •  3 

1-5 

1-7 

£  -9 

-1.02772(-1) 

-8.54582(-2) 

5 

-A.18326(-2) 

-3.13956(-2) 

6 

-2.48972(-2) 

-2.00775(-2) 

-1.38983(-2) 

7 

-1.64735(-2) 

-1.39314(-2) 

-1.02101(-2) 

8 

-1.16248(-2) 

-1.02309(-2) 

-7.81653(-3) 

-7.81419(-3) 

9 

-8.53928(-3) 

-7.83020(-3) 

-6.17575(-3) 

-6.17407(-3) 

10 

-6,65276(-3) 

-6.18573(-3) 

-5.00213(-3) 

-5.00092(-3) 

-5.00056(-3) 

11 

-5.32418(-3) 

-5.009l7(-3) 

-4.13387(-3) 

-4.13296(-3) 

-4.13267(-3) 
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Table  VIII.  Energy  levels  of  type  Rydberg  states*. 


n 

f  -  2 

1-4 

1-6 

00 

N 

U 

-5.10714(-2) 

5 

-2.85280(-2) 

-2.00253(-2) 

6 

-1.83158(-2) 

-1.39036(-2) 

7 

-1.27814(-2) 

-1.02134(-2) 

-1.02075 (-2) 

8 

-9.43470(-3) 

-7.81880(-3) 

-7.8l484(-3) 

9 

-7.25213(-3) 

-6.17733(-3) 

-6.l7453(-3) 

-6.17373(-3) 

10 

-5.74681(-3) 

-5.00330(-3) 

-5.00126(-3) 

-5.00066(-3) 

11 

-4.65981(-3) 

-4.13475(-3) 

-4.l3322(-3) 

-4.13277(-3) 

a 


See  the  footnotes  in  Table  VI 
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1  a 

Table  IX.  Energy  levela  of  type  Rydberg  states  • 


n 

t  -  1 

1-3 

1-5 

1-7 

1-9 

3 

-9.13323(-2) 

A 

-A.408A6(-2) 

-3.13501(-2) 

5 

-2.60713(-2) 

-2.00529(-2) 

6 

-1.72373(-2) 

-1.39201(-2) 

-1.38974(-2) 

7 

-1.22AA2(-2) 

-1.022A0(-2) 

-1.02095(-2) 

8 

-9.lA621(-3) 

-7.82596(-3) 

-7.816l5(-3) 

-7.81411(-3) 

9 

-7.09158(-3) 

-6.18238(-3) 

-6.17547(-3) 

-6.17401(-3) 

10 

-5.65876(-3) 

-5.00700(-3) 

-5.00193(-3) 

-5.00090(-3) 

-5.00054(-3) 

11 

-A.61911(-3) 

-4.1375A(-3) 

-4.13373(-3) 

-4.13293(-3) 

-4.13266(-3) 

^  See  the  footnotes  In  Table  VI. 


3  +  a 

Table  X.  Energy  levels  of  E  type  Rydberg  states  . 

E 


n  £-2  £-6  £  =  8 


3 

4 

5 

6 

7 

8 
9 

10 


-1.13895(-1) 

-4.86088(-2) 

-2.73442(-2) 

-1.75921(-2) 

-1.229l4(-2) 

-9.08172(-3) 

-6.98621(-3) 

-5.53839(-3) 


-5.79860(-2) 

-3.23170(-2) 

-2,05351(-2) 

-1.41883(-2) 

-1.03858(-2) 

-7.92941(-3) 

-6.25120(-3) 

-5.05343(-3) 


-2,00298(-2) 

-1.39063(-2) 

-1.02151(-2) 

-7.81988(-3) 

-6.71812(-3) 

-5.00386(-3) 

-4.13518(-3) 


-l,02078(-2) 

-7.81499(-3) 

-6.17465(-3) 

-5.00134(-3) 

-4.13328(-3) 


-6.17376(-3) 

-5.00067(-3) 


4.48675(-3) 


-4.13278(-3) 
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3  4>  A 

Table  XI.  Energy  levels  of  type  Rydberg  states  . 


t  »  1 

1-3 

1-5 

1-7 

l.03570(-l) 

4,25219(-2) 

2.52225(-2) 

-8.72696(-2) 

-3.14252(-2) 

-2.00996(-2) 

-1.38983(-2) 

1.66328(-2) 

-1.39399(-2) 

-1.02101(-2) 

1.17026(-2) 

-1.02361(-2) 

-7.81654(-3) 

-7.81420(-3) 

■8.87308(-3) 

-7. 832 64 (-3) 

-6.17575(-3) 

-6.17407(-3) 

•6.66733(-3) 

-6.18715(-3) 

-5.00214(-3) 

-5. 00093 (  3) 

•5.33223(-3) 

-5,01050(-3) 

-4.13388(-3) 

-4.13297(-3) 

-5. 00056 (-3) 
-4.13267(-3) 


y  ^  U  V  ^ 

^•Vj 


4 

-5.16214(-2) 

5 

-2.87212(-2) 

-2.00254(-2) 

6 

-1.84182(-2) 

-1.39037(-2) 

7 

-1.28411(-2) 

-1.02135(-2) 

-1.02075(-2) 

8 

-9.46553(-3) 

-7.81880(-3) 

-7. 81484 (-3) 

9 

-7.27266(-3) 

-6.17733(-3) 

-6.17453(-3) 

-6.17371(-3) 

10 

-5.76093(-3) 

-5.00331(-3) 

-5.00126(-3) 

-5.00066(-3) 

11 

-4.66932(-3) 

-4.13475(-3) 

-4.13322(-3) 

-4.13277(-3) 

®See  the  footnotes  in  Table  VI. 
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Table  XIII. 


Energy  levels  of 


type  Rydberg  states*. 


n 

i  »  1 

1  -3 

t  -  5 

1-7 

1-9 

0 


3 

-9.19470(-2) 

U 

-A.A2717(-2) 

-3.13516(-2) 

5 

-2.61A97(-2) 

-2.00579(-2) 

6 

-1.72766(-2) 

-1.39237(-2) 

-1.38974 (-2) 

7 

-l,22662(-2) 

-1.02265(-2) 

-1.02095(-2) 

8 

-9.15955(-3) 

-7.82781(-3) 

-7.81615(-3) 

-7.8l4ll(-3) 

9 

-7.10009(-3) 

-6.18374(-3) 

-6.17548(-3) 

-6.17401(-3) 

10 

-5.66435(-3) 

-5.00803(-3) 

-5.00193(-3) 

-5.00090(-3) 

-5.00054(-3) 

11 

-A,62274(-3) 

-4.13833(-3) 

-4.13373(-3) 

-4.13293(-3) 

-4.13266(-3) 

a 


See  the  footnotes  in  Table  VI. 


,%  .V 

Sv*** 


m 


>  r* 


I 


y  ^ 

n  =  4 

-0.749 

0.624 

n  =  5 

-0.216 

-0.545 

0.770 

0 

n  =  6 

-0.188 

-0.435 

0.837 

0 

n  =  7 

0.159 

0.368 

-0.875 

-0 

n  =  8 

-0.140 

-0.321 

0.902 

0 

n  =  9 

-0.124 

-0.281 

0.924 

0 

n  ^  10 

-0.109 

-0.245 

0.944 

0 

n  =  11 

0.092 

0.202 

-0.970 

a 


Jl  =  2  in  all  basis  functions  (k,0 


^  -  fc 
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Table  XIX,  Comparison  of  energy  levels  with  other  works- 


Present 


Ref.  3‘ 


Experiment 


12.785 

13.1 

12.98 

12.935 

1A.4A3 

14.4 

14.35 

14.364 

14.904 

14.9 

14.86 

14.871 

15.133 

15.13 

15.110 

15.265 

15.245 

15,349 

15.332 

15.400 

15 . 386 

15.436 

15.426 

13.096 

13.1 

12.95 

12.982 

14.381 

14.5 

14.32 

14.329 

14.872 

15.0 

14.81 

14.846 

15.112 

15.11 

15 .098 

15.248 

15.238 

15.332 

15.327 

15.388 

15.385 

15.427 

15.425 

15.455 

15.454 

12.689 

12.5 

12.25® 

14.318 

14.0 

14.191 

13.06 

14.100 

12.763 

13.0 

12.84^ 

14.424 

14.4 

14.895 

14,9 

13.079 

13.0 

12.8® 

14.376 

14.4 

14.869 

14.9 

12.482 

11.9 

11.87® 

14.258 

13.9 

Measured  in  eV  above  N  (X  E  ) ;  See  text. 

2  g 

^Singlet-triplet  distinction  cannot  be  made;  see  text 


^Ref.  17. 
g 

®Ref,  1. 


Ref. 18. 


PART  II 

OPTICAL  EXCITATION  CROSS  SECTION  FOR  EMISSION  OF 
LONG-WAVE  INFRARED  LINES  OF  ATOMIC  OXYGEN  PRODUCED 
BY  ELECTRON  IMPACT  ON  OXYGEN  MOLECULES 


Long-wave  infrared  emission  of  atmospheric  atoms,  especially  oxygen,  is  of 
great  interest.  One  important  mechanism  for  producing  the  "oxygen  long-wave 
infrared  emitters",  l.e.,  highly  excited  oxygen  atoms,  is  dissociation  of  O2  by 
electron  Impact, 

O2  0  ♦  0». 

The  ground  state  oxygen  atom  has  an  electron  configuration  Is  2s  2p  .  The 
excited  states  which  lie  below  the  first  ionization  limit  are  well  described  by 
the  Russell-Saunders  coupling.  A  configuration  2p^nt  comprises  a  number  of 
levels  characterized  by  L,  S,  and  J.  In  this  work  we  are  not  concerned  with  the 
fine  structure,  thus  we  designate  an  excited  level  by  n,  t,  L,  S. 

An  oxygen  atom  in  a  highly  excited  level,  say  n=7,  may  emit  long-wave 
infrared  radiation  upon  decaying  to  a  level  of  slightly  lower  energy,  say  n=6. 

For  example  the  7f  +  6d  transition  gives  emission  at  120,577  A.  However, 
measurement  of  the  intensity  of  this  emission  is  very  difficult  because  the 
infrared  detectors  available  for  120,000  A  wavelength  have  very  low  sensitivity 
compared  to  photomultipliers  which  operate  in  the  wavelength  range  of  about 
2,000  -  12,000  A. 

Recently  Schulman  et  al.^  have  performed  experiments  on  production  of 
excited  oxygen  atoms  by  electron  impact  on  oxygen  molecules  with  incident 
electron  energy  from  threshold  to  500  eV.  They  measured  the  intensity  of 
radiation  from  the  excited  oxygen  atoms  in  the  wavelength  range  3,690  -  11,300  A 
corresponding  to  transitions  from  the  highly  excited  states  to  the  low  excited 
states  (3s  and  3p).  Absolute  optical  emission  cross  sections  have  been  reported 
for  transitions  in  the  wavelength  range  3,690  -  11,300  A  originating  from  some 
thirty  spectroscopic  terms.  In  this  work  we  show  that  by  combining  the  data  of 
Schulman  et  al .  with  theoretical  transition  probabilities,  it  is  possible  to 
obtain  optical  cross  sections  for  emission  of  atomic-oxygen  lines  in  the 
long-wave  infrared  region  (produced  by  electron  impact  on  oxygon  molecules)  which 
arc-  otherwise  very  difficult  to  measure  directly. 
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The  ratio  of  the  optical  ealsslcn  croaa  sections  for  two  transitions  which 
originate  from  the  same  upper  level  is  equal  to  the  ratio  of  the  corresponding 
transition  probabilities  A,  l.e. , 


(1) 


Thus  optical  cross  sections  for  the  emission  lines  that  share  a  common  upper 
level  with  a  transition  reported  in  Ref.l  can  be  determined  by  using  the 
appropriate  transition  probabilities.  Since  direct  measurement  of  emission 
intensity  in  the  Infrared,  especially  long-wave  infrared,  region  Is  more 
difficult  Chan  in  the  visible  and  near  ultraviolet  region,  Eq.  (l)  allows  us 
to  obtain  optical  cross  sections  for  infrared  lines  from  the  theoretical 
transition  probabilities. 

In  this  section  we  consider  only  transitions  between  levels  with  the 

3  4  +4 

2p  (  S)nl.  configuration.  Since  the  0  core  is  in  a  S  staCr,  the  total  orbital 

angular  momentum  L  of  the  oxygen  atom  is  equal  to  t.  An  excited  level  is 

then  designated  by  nLSJ.  The  transition  probability  between  two  levels 
2 

may  be  written  as 


A(n'LSJ'  -*■  nLSJ)  (M)j^(L)a(n'L'S  nLS)/(2J’+l) , 

a(n*L’S  -►  nLS)  »  64ir^'^^a^/3hc^, 

-  e^|/R^j^(r)rR^,j^,(r)r^dr|*/(4l^-l), 


(2) 

(3) 

(4) 


wherejfClO  is  the  multiplet  strength, kf  (L)  the  line  strength,  R  .  the  radial 

nJc 

wave  function,  and  is  the  greater  of  i  and  V.  Formulas  forJii(H)  and 

J(L)  are  given  in  Ref.  3.  For  transitions  between  states  of  the  type 

2p^(^S)ni,  the  values  ofj«|(M)  are  15  for  ^S-^P,  150  for  ^P-^D,  525  for  ^D-^F, 
3  3  3  3  3  3 

9  for  S-  P,  90  for  P-  D,  and  315  for  D-  P.  The  line  strengths  are  shown 
in  Table  1. 
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The  radial  wave  functions  R^(r)  in  Eq.  (A)  are  calculated  by  the 

Hartree-Fock  self-consistent-field  (SCF)  procedure.  To  sinplify  the 

calculation,  the  Is  and  2s  orbitals  of  the  excited  oxygen  atom  are  taken 

+  3  4 

from  the  SCF  solution  of  the  0  [2p  (  S)]  ion,  and  are  not  varied  in  the 

Iteration  cycles.  In  Table  II  we  list  the  values  of  a(n'L'S  nLS)  calculated 

from  these  Hartree-Fock  wave  functions  for  the  transitions  which  are  used 

to  obtain  the  long-wave  Infrared  emission  cross  sections  according  to  Eq.(l). 

3  4 

For  a  given  configuration  2p  (  S)nt,  the  radial  wave  functions  for  the  triplet 
states  are  different  from  those  for  the  quintet  states  because  of  the  different 
exchange  terms  in  the  Fock  equations.  This  Is  particularly  true  for  the  np 
series  on  account  of  the  strong  exchange  Interaction  between  the  2p  and  np 
orbitals.  We  see  in  Table  u  that  this  difference  in  radial  functions  between 
the  triplet  and  quintet  series  Is  reflected  In  the  transition  probabilities. 

To  get  an  indication  of  the  accuracy  of  the  Hartree-Fock  calculation, 
we  have  computed  the  transition  probabilities  using  a  different  method,  i.e., 
the  self-lnteract Ion-corrected  form  of  the  local-spln-denslty  (LSD)  approximation. 
This  method  has  been  discussed  extensively  in  the  llterature,^~^  thus 
only  a  brief  description  is  presented  here.  In  the  conventional  LSD  or  the 
Hartree-Fock-Slater  method,  the  exchange  potential  in  the  Fock  equations  for 
the  one-electron  orbitals  is  approximated  by  a  local  potential  proportional 
to  the  cubic  root  of  the  electron  density  (on  spin  density  for  spin-polarized 
cases),  p  .  This  local  exchange  approximation  greatly  simplifies  the 
computational  work,  but  has  a  serious  dranbaek  related  to  the  self-interaction 
in  the  following  way.  The  effective  potential  of  an  electron  in  the 
Hartree-Fock  theory  consists  of  the  electron-nucleus  attraction,  electron- 
electron  Coulomb  repulsion,  and  the  exchange  potential.  The  electron- 
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electron  Coulomb  repulsion  Includes  a  term  corresponding  to  the  interaction 

of  the  electron  with  itself.  This  unphysical  self-interaction  is  cancelled 

by  an  identical  term  in  the  exchange  potential  so  that  its  presence  has  no 

effect.  However,  when  the  exchange  potential  is  approximated  by  the  local 
1/3 

p  form,  the  cancellation  of  the  self-interaction  energy  is  not  complete. 

One  troublesome  consequence  of  this  incomplete  cancellation  is  that  for  a 

neutral  atom  the  electron  potential  at  a  large  distance  r  from  the  nucleus 

approaches  zero  much  faster  than  the  correct  -1/r  dependence,  resulting  in 

an  underestimation  of  the  net  attraction,  hence  the  ionization  energy. 

This  difficulty  Is  partly  resolved  by  a  cutoff  procedure,  adopted  in  the 

8 

Herman-Sklllman  formulation  of  the  Hartree-Fock-Slater  method,  in  which  the 

local  exchange  approximation  is  used  only  in  the  region  of  r  between  zero 

and  a  cut-off  value  r  ,  and  for  r>r  the  electron  potential  is  simply  set  to 

o  o 

-1/r.  Nevertheless  this  cutoff  procedure  does  not  completely  remove  the 
spurious  self- Interact Ion,  although  it  gives  the  correct  asyrptotlc  potential 
energy. 

A  more  fundamental  way  of  addressing  the  problem  of  self-interaction  is 

to  remove  the  self -interaction  terms  from  the  total  energy  in  the  outset. 

This  approach  has  recently  met  with  great  success  in  calculations  of 

4-6 

electronic  structure  of  atoms,  and  is  often  referred  to  as  the  self- 

interaction  correction  (SIC).  Of  special  interest  to  us  is  the  paper 
of  Harrison  et  which  demonstrates  the  success  of  the  SIC-LSD  approximation 

for  calculating  energies  of  excited  states  of  atoms.  The  SIC-1.SD  approximation 
retains  the  computational  simplicity  of  the  Hartree-Fock-Slater  approach, 
but  provides  a  significant  improvement  over  the  conventional  l.SD  approximation. 


9 


and  the  results  agree  well  ('''10%) 


with  the  values  obtained  by  the  Hartrce-Fock  wave  functions  (Table  ii). 

The  agreement  between  the  two  sets  of  results  supports  the  accuracy  of  eur 

calculation  within  the  Hartree-Fock  framework. 

Earlier  works  on  the  transition  probabilities  of  the  oxygen  atom 

include  the  calculation  using  wave  functions  determined  by  the  Hartree-Fock- 

Slater  method  (without  the  SIC).^^  The  formulation  of  the  Hartree-Fock-Slater 

method  as  given  in  Ref.  8  uses  the  total  electron  density  (rather  than  the 

1/3 

spin  density)  In  the  p  -exchange  term;  therefore  it  gives  the  same 

3  4 

R  o(r)  for  the  quintet  and  triplet  states  of  the  2p  (  S)nS,  configurations, 
nt 

Nevertheless,  compared  with  our  values  of  O  [Eq.  (4)  1  of  the  quintet  series, 
those  of  Ref. 10  are  generally  within  a  25Z  range  when  Q  Is  greater'than  1.0. 
For  smaller  values  of  a  the  discrepancy  is  much  larger.  As  pointed  out  In 
Ref.  10,  the  integral  of  Eq.  (4)  Is  very  sensitive  to  the  wave  functions, 
particularly  when  severe  cancellation  occurs  between  the  positive  and 
negative  contributions.  Thus  it  is  not  surprising  to  find  a  much  larger 
discrepancy  between  the  Hartree-Fock  and  Hartree-Fock-Slater  results  for  the 
cases  where  Q  Is  small. 

Blemont  and  Grevesse^^  calculated  oscillator  strengths  for  a  large 
number  of  Infrared  lines  of  atomic  oxygen  using  the  Coulomb  approximation. 

The  Coulomb  approximation  takes  advantage  of  the  fact  that,  in  the  usual 
length  formula,  nearly  all  the  contribution  to  the  dipole  matrix  element  comes 
from  the  outer  lobes  of  the  two  wave  functions  Involved.  This  asymptotic 
part  of  the  wave  function  may  be  computed  simply  with  the  knowledge  of  the 
Coulombic  potential  and  the  ionization  energy,  since  the  short-range 
potentials  such  as  electron-exchange  potential  become  negligible  for  large 
r.  Therefore,  this  method  Is  suitable  for  computing  the  dipole  matrix 
elements  between  a  pair  of  highly  excited  states.  It  follows,  on  the  other 
hand,  that  the  validity  of  this  method  may  be  questioned  for  low-lying 


51 

states,  and  this  point  is  discussed  in  Ref.  12. 

The  majority  of  the  infrared  lines  of  Ref. 11  involve  highly  excited 

2 

states.  Comparison  is  made  of  O  in  Eq. (4)  between  the  present  work  and 

Ref .11  for  some  140  transitions  (ns,  np,  nd,  and  nf,  with  n  up  to  10). 

In  most  cases  the  agreement  is  within  25Z.  About  ten  cases,  in  which 

3  3  3  3 

discrepancy  Is  greater  than  25Z,  occur  mostly  in  S  -  P  and  P  ~  D  series. 

In  computing  dipole  matrix  elements,  the  relative  "phase"  of  the  pair 
of  wave  functions  is  an  important  factor.  The  position  of  the  outer  peak 
of  a  wave  function  is  dictated  by  the  details  of  the  potential  in  the 
interior  region.  In  our  Hartree-Fock  calculation  the  triplet  and  quintet  func¬ 
tions  of  the  same  nf.  are  shifted  relative  to  each  other  because  of  the 
difference  in  exchange-coupling  coefficients.  When  the  Coulomb  approximation 
is  used,  such  effect  is  taken  into  consideration  to  some  extent  through  the 
effective  quantum  number  n*,  but  falls  short  of  the  rigorous  treatment  with 
the  result  that  the  wave  function  may  be  shifted  by  an  improper  amount  from 
the  hydrogenlc  counterpart.  Although  the  inaccuracy  in  the  phase  shift  may 
be  small  in  each  of  the  wave  functions,  it  could  lead  to  much  greater  error 

in  the  dipole  matrix  elements  where  the  relative  phase  shift  of  the  two 

3  3  3  3 

functions  is  an  Important  factor.  The  discrepancy  in  the  S  -  P  and  P  -  D 
series  referred  to  in  the  preceding  paragraph  may  be  attributed  to  the 
inaccuracy  in  the  phase  shift  in  the  Coulomb-approximation  wave  functions. 

The  values  of  a(n’L’S  nLS)  in  Table  n  are  combined  with  the  measured 
optical  emission  cross  sections  of  Ref.l  in  accordance  with  Eq.  (l)  to 
yield  optical  cross  sections  for  a  number  of  infrared  lines.  Since  the 
individual  J'  -►  J  components  of  an  n'L'S  nLS  transition  are  not  resolved 
in  the  measured  emission,  we  assume  that  the  population  of  the  J-levels  (due 
to  electron  impact)  within  an  (nLS)  term  is  proportional  to  the  statistical 
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Table  II.  a(n*L*S  ^  nLS)  as  defined  In  Eq.  (3). 


0  MulClplet 

X(to 

a  V  sec  ) 

6s  -►  3p 

5437. A5 

1.28824  (6)  * 

3  3 

6s  S  3p  P 

6048.06 

1.05174  (6) 

^  5^  ,  5„ 

6s  S  -►  5p  P 

71779.87 

5.44747  (5) 

3  3 

6s  S  ^  5p  P 

77211.17 

4.49093  (5) 

7s  h  ^  3p  ^P 

5020.83 

7.12077  (5) 

7s  ->  6p  ^P 

132148.12 

2.09108  (5) 

4p  -»■  3s 

3948.57 

3.50477  (5) 

4p  3d 

59773.34 

4.93382  (4) 

5d  3p  ^P 

5331.32 

4.54278  (5) 

5d  4f 

41371.10 

2.85093  (3) 

5d  -►  5p  ^P 

56842.49 

2.76437  (5) 

7d  ->  3p  ^P 

4774.37 

1.17938  (5) 

7d  -»■  5f 

46907.72 

1.09355  (3) 

7d  -►  6p  ^P 

56026.60 

4.03733  (4) 

7d  -►  6f  ^F 

126850.42 

1.39241  (3) 

6f  ^F  ->  3d 

10678.77 

1.57035  (5) 

3  3 

6f  F  3d  D 

10756.48 

1.50984  (5) 

6f  ^F  ^  5d 

71440.57 

5.01200  (4) 

3  3 

6£  F  5d  0 

72664.23 

5.12763  (4) 

7f  ^F  3d 

9828.63 

8.86534  (4) 

7£  ^F  3d  \ 

9894.46 

8.47534  (4) 

7f  ^F  5d 

45254.09 

3.01991  (4) 

3  3 

7£  F  -*•  5d  D 

45742.75 

3.01714  (4) 

7£  ^F  >  6d 

118529.35 

1.76778  (4) 

3  3 

7f  F  »■  6d  D 

120577.23 

1.83095  (4) 

a 

Kumbers  in  the  parentheses 

denote  the  power 

o£  10. 

Table  HI*  Excitation  cross  sections  Cor  long-wave  infrared  emission 
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of  0  produced  by  dissociative  electron- impact  excitation  of  O2  at  100  eV. 


0  IKiltiplet 


x(X) 


5d 

a 

0 

4f 

7f 

5d 

^0° 

7f 

-► 

5d 

^0° 

7d 

-► 

5f 

7d 

-V 

6p 

5d 

5p 

4p 

-► 

3d 

^0° 

6f 

5d 

6s 

-*■ 

5p 

6f 

^F 

5d 

6s 

-► 

5p 

7f 

6d 

^0° 

7f 

3 

F 

-► 

6d 

7d 

6f 

^F 

-► 

6p 

41371 

45254 

45743 

46908 

56027 

56842 

59773 

71441 

71780 

72664 

77211 

118529 

120577 

126850 


(10"^°  cm^) 


0.075 

0,82 

0.53 

0.019 

0.20 


2.1 


2.9 


1.4 


0.33 

1.3 


0.13 

0.48 

0.32 

0.024 


132148 


0.085 
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